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The ma 6 (m/M) order corrections to the hyperfine splitting in the ion are calculated. That 
allows to reduce uncertainty in the frequency intervals between hyperfine sublevels of a given rovi- 
brational state to about 10 ppm. Results are in good agreement with the high precision experiment 
carried out by Jefferts in 1969. 



I. INTRODUCTION 

In our previous work [l[ we have calculated the hyperfine structure of the hydrogen molecular ion Hj within the 
Breit-Pauli approximation taking account of the anomalous magnetic moment of an electron. This approximation 
includes the contributions of order ma 4 (m/M) and ma 5 (m/M) and thus the relative uncertainty in determination 
of the hyperfine structure intervals is of about 5 x 10~ 5 . For the first time that has allowed to confirm the Jefferts 
measurements [2| to the level of experimental accuracy of 1.5 kHz for transitions within the same multiplet F (F is 
the total spin of a state in ion). For the spin-flip transitions (F = 3/2) — » (F = 1/2) a discrepancy of about 80 kHz 
still remains. 

The main goal of the present work is to consider higher order corrections to the hyperfine splitting of H^ to reduce 
the discrepancy with the Jefferts experiment for spin-flip lines down to a few ppm. To that end we will calculate the 
QED contributions of order a 2 Ep and partially of order o?Ep along with the proton finite size corrections such as 
Zcmach and pure recoil contributions, which are essential at this level of accuracy. 

The effective Hamiltonian of the spin interaction for the H^ ion is (we use notation of [l|): 



H cS = b F (I ■ s e ) + c e (L • s e ) + c/(L • I) 



di 



(2L-l)(2L+3) [3 
d 2 



(2L-l)(2L+3) 



L 2 (I.s e )-[(L.I)(L-s e ) + (L-s e )(L-I) 



1l 2 I 2 -1(L.I)-(L.I) 2 



(1) 



here I is the total nuclear spin, L is the total orbital momentum. The assumed coupling scheme of angular momenta 
is: P = I + s e , J = L + P. 

The major coupling is the spin-spin electron-proton interaction (first term in {1}) which determines the principal 
splitting between F — 1/2 and F = 3/2 states. So, the main contribution to the theoretical uncertainty on the 
spin- flip transition frequencies is uncertainty in the spin-spin interaction coefficient b p , and our aim is to calculate an 
improved value for bp. 

Here it is very useful to make a comparison with the HFS studies of the hydrogen atom ground state. Indeed, the 
analytical form of many contributions to the hyperfine splitting of H^ can be obtained from these results. Moreover, 
we will use the known results on the hydrogen atom as a guide and a check of our analytical derivations. 

The hyperfine splitting for the ground state of a hydrogenlike atom may be obtained with high accuracy already 
from the nonrelativistic quantum mechanics (see for example Q), 



16 2 m e 
E F = -ja cR^ p — 



-3 



(2) 



here fi p is the magnetic moment of a proton in nuclear magnetons, m e and M p are the electron and proton masses, 
respectively. Quantum electrodynamics corrections without recoil terms have been known for some time 0, [E| and 
may be expressed as: 



hfs,QED 



E F 



■(Za) 



In 2 



a(Za) 



3tt 



a{Za) 2 \n 2 {Za) 



(3) 
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term 


/ITT N 

(kHz) 


E F 


1418 840.09 


a e E F 


1645.361 




113.333 




-136.517 


AEat(Za) 2 ln 2 (Za) 


-11.330 


higher order QED 


1.26 


AE Z 


-56.9(9) 


AE* R 


8.43(8) 




2.0(8) 


A£(HFS) 


1420 405.7(1.7) 


experiment 


1420 405.7517667(9) 



TABLE I: Contributions to the hyperfine splitting of the ground state in a hydrogen atom. Uncertainty in the Fermi splitting, 
Ef, is determined by the fine structure constant. The second row is the contribution due to the anomalous magnetic moment 
of an electron, a e . Contributions of higher orders, (Za) 2 , a(Za), etc. as well as proton finite size corrections, are taken from 
Eqs. 



where a e is the electron anomalous magnetic moment. We keep Z, the nuclear charge number, in all expressions in 
order to make clear the origins of different corrections. 

Beyond pure QED corrections there are also proton structure effects (see [1, [H, 0] for a detailed discussion). The 
leading one is the Zemach correction [H ((Za)(m/A)Ep) that along with the radiative corrections to the nuclear 
structure contribution [9(] reads 

Tfl A/f 

AE Z = -2 e p (Za) R Z (1 + 5% d )E F , (4) 
Mr, + m P 



Here Rz is the Zemach radius, a mean radius associated with the proton's charge-current distribution, 
1 f d 3 q 

7 



1 G E (-q 2 )G M (-q 2 



where Ge and Gm are the electric and magnetic form factors of a proton. We take Rz = 1.045(16) fm Q. Radiative 
corrections 8^ d to the Zemach contribution were obtained in [9( and (5™ d = 0.0153. The parameter A determines the 
energy scale that corresponds to the mean radius of the proton, and A « 0.8Af p Q. Next are the pure recoil proton 
structure corrections of orders (Z a) k (m/ M)Ef (k = 1,2) [l(| 

ks (5.84 ± 0.01) x 10 -6 E F . (5) 

The last remaining effect of order (Za)(m/M)Ep , which has to be included, is the proton polarizability [TTI ] 

A pol w (1.4±0.6)xlO" 6 £; F . (6) 

A summary of various contributions to the HFS of the hydrogen ground state is given in Table H Up to now, m our 
previous studies for the H J ion [l[ , only the contributions from the first two lines have been included into consideration. 
In the present work we intend to extend our research to higher order QED corrections (up to Ep a(Z ot) 2 \x? [Z a) 
term) as well as the proton structure effects. 

The major part of the contributions mentioned above in Eq. ((4j- ([6j) may be considered as contact type interactions, 
which depend on the value of the squared density of the nonrelativistic wave function at the electron-proton coalescence 
point. Thus they do not require new extensive calculations, the mean values for the delta function operators can be 
taken from The main task is calculation of the relativistic correction term of order Ep(Za) 2 , which may 

be performed using the nonrecoil limit of the two center problem. The obtained effective adiabatic potentials are 
subsequently averaged over the radial wave function as it was done for the ma 6 order relativistic correction to ro- 
vibrational energies in IB 

The paper is organized as follows. In Sec. II. A and II. B we use the NRQED to derive all the spin-dependent 
interactions of order a 6 (m/M) and the corresponding potentials in the coordinate space. The radiative corrections 
of orders a(Za)Ep and a(Za) 2 hi 2 (Za)Ep, as well as proton structure effects, which may be expressed as contact 
type interactions in the NRQED, are given in paragraph C. The perturbation formalism used to obtain the energy 
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corrections is described in paragraph D. In the next two sections, we describe the calculation of the relativistic 
corrections of order (Za) 2 Ep. First, for the HFS of the hydrogen ground state (Sec. Ill), where we rederive the well- 
known Breit correction [15| , providing a useful check of our approach. Then in Sec. IV the H J ion case is considered. 
Finally, numerical results are given and discussed in Sec. V. 



II. NRQED INTERACTIONS 



In this section we use the NRQED [17[ to describe the interactions, which are of relevance to our problem. A nice 
and illuminative introduction to the NRQED approach may be found in [f|. The units c = % = 1 and e 2 — a are 
used in this section, the elementary charge, e, is positive. We consider the low energy scattering, assuming that the 
momentum of a particle is of order Za, and we expand the scattering amplitude in terms of a and p 2 . 



A. Tree-level interactions of order ma 6 (m/M). 
The momentum 4- vectors for the scattering of an electron (or proton) by the field of a static external source obey: 



Po=Po = E i 



qo = 0, 



q 



where p and p' are 4-moments of incident and scattered particles, respectively. 

On-shcll Dirac spinors can be presented via the Schrodingcr-Pauli spinors as follows 



u(p) 



I E„ 



X 



2E V 



Z p 

, E p +m 



X 



u*{p)u{p) = (X*X) = 1, 



E p = y/m 2 



p 2 = m I 1 



2m 2 8m 4 16m 6 



here erf are the two component Pauli matrices and X are the two-component Schrodinger-Pauli wave functions. We 
assume that Dirac spinors are normalized as (u*u) = 1. That corresponds to the nonrelativistic normalization: the 
probability to discover a particle in a unit of volume is equal to unity. With this normalization Dirac spinors are 
expanded in the low-energy limit as follows 



u{p) 



1 



p 



Z P 

2m 



Hp" 

128m 4 

3pI 
sni- 



X' 



X 



The nonrelativistic scattering amplitude at tree-level for a scalar static field is determined by the following expansion 



A E (p,p')=X*(p')(eZ)A (d) 1 



. <r p [q x p] 3q 2 (p' : 



8m 2 



4m 2 



64m 4 



128m 4 

and for a vector static field one obtains 



Hp' 2 -P 2 ) 2 _ .3<r p [q x p](p' 2 + p 2 ) 



32m 4 



(7) 



X(p) 



A ( '\ Y*f n '\(~7\At n \l P , + P .Kxq] , P'(3y 2 +p 2 )+ P (3p 2 +y 2 ) 

Am{P,P ) = X (p )(eZ)A(q) i- 



+i 



2m 2m 16m 3 

[<t p x p'](3p' 2 + p 2 ) - [a p x p](3p 2 + p' 2 ) 
16m 3 



(8) 



where Z is the charge of a particle e, p,. . . For an electron Z = — 1. 

In what follows an index a or b = 1,2 denotes nucleus 1 or 2 in the ion, indices i,j = 1,3 are Cartesian 
coordinates. The imaginary unit is denoted by upright i. 

The higher order vertices of tree level diagrams produce new interactions (q = p' e p e ): 
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a) via Coulomb photon exchange: 

.3crf[qxp e ](y e 2 +p2)^ 1 
|/i = e 1 



32mf 

b) via transverse photon exchange: 

. [«rf x p'J (3^ 2 +p 2 ) - [a* x p e ] (3p 2 +^ 2 ) ^ 



V 2 = e z -i 



V 3 = e A 



p;(3p?+ p 2) + Pe ( 3 p2 + J , ?r 



-Xts 



— — 5 lJ 



q l q> 



-Z a 



Pq + Pq 

2M a 



V 4 = e 2 -i 



16m^ 



>f xp' e ](3p' e 2 +p 2 e ) - [of xp e ](3p 2 +?a 



■ 7 of x (-q) 



2M a 



5« 



• 7 of X (~q) 

2M a 



(9) 

(10a) 
(10b) 
(10c) 



In parentheses here are the vertex functions of the effective NRQED interaction taken from Eqs. (0)-® • The approx- 
imate transverse photon propagator (see [Hj], § 83) is placed in the square brackets. 
The obtained potentials can be simplified as follows: 



V 2 = e 2 -i 



2 / ; 2[of xq](p' e 2 +p 2 ) + [of X (p' e + p e )](p' e 2 -p 2 ) 



16m;? 



The last term in the first brackets produces a symmetric operator with the property: {(p, Aip) = 0, for an arbitrary if. 
That means that this operator is identical to zero operator, and V2 may be rewritten 



V, = -e 2 i 



[<rf xq](^ 2 +p 2 )V 1 



8m2 



~ {Za w n 



In a similar way the other operators may be simplified: 



V, = e 2 



V 4 



Pejpf+PP Y 1 ( ofxq 
4m| / q 2 V ° 2M « 



[of xq](p' e 2 +p 2 ) V 1 of xq 



8m5! 



2M„ 



(11a) 



(lib) 



(11c) 



Transforming potentials V„ to the coordinate space (r a = R e — R a , where R e , R a are the coordinates of electron 
and nuclei with respect to the center of mass) one gets, using the notation {A, B} — AB + BA: 



V, = -a 



■iz a 



16rMp I r; 



Pl> -4[raXPe] \ S 



Vo = a 



Z a 



tm\M a 



Pe>i[r«xP a ] \s e , 



1 



V3 = - a ^3\ P *>^ raXPe ^^> 



(12) 



V 4 = -a 



4m;? 



p 2 , 



8?r „ v _ rgs e /A a -3(s e r a )(/x a r a ) 
s e ^z a 0(r a< ) 



Here /i a is the magnetic moment operator for nucleus a. Only V4 involves both electron and nuclear spins and 
contributes to bp- 
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o 



/ X 



X X 

e 2 E 2 





8m 3 

6 



e 2 o- ■ [A x E] 

4m 2 

.6 



e 2 A 2 
2m 



FIG. 1: Top row: QED scattering diagrams. Bottom row: NRQED seagull interactions inferred from QED diagrams. 



B. Seagull-type interactions 



In Fig. [T] three NRQED seagull diagrams are presented. They may be obtained from the corresponding QED 
Z-diagrams by expanding the scattering amplitude in terms of p 2 . The double Coulomb photon exchange diagram 
has a leading order a 6 , however it does not involve interactions dependent on spin. The third diagram is double 
transverse photon exchange has a recoil order (m/M) 2 , and is also out of interest for present consideration. 

The potentials which stems from the seagull vertex with one Coulomb and one transverse photon lines can be 
expressed as follows (qi = Pi P' l5 q2 = P2 P^) 1 



V, 



V 6 




4 



-z h 



-iZ, 



Pi 



2M b 



x 1 



qo 



x (-qt) 

2M b 



q« 
q« 



The sources Z a and ifb may belong to a same particle or to two different particles. 
In the coordinate space one has, for a ^= 6, 



V s = -a 



Z a Z b f [r x P fc ] [r x r b ] (r 6 P 6 ) 



^6 



V 6 



4t7J 2 : M 6 

2 Z a ZbHb [r a x s e ] [r b x s b ] 
2m?M 6 ' 



r 3 r 3 



a 



Z a [r a x s e ] [r 6 x /x b ] 



2mi 



r 3 r 3 



Using [axb][cxd] = (ac)(bd) — (cb)(ad), one may further simplify 
When the sources coincide (a = b), the interactions modify as follows: 



V 7 = -a' 



Z 2 a [r n xP g ] : 

AmlMa r a 

Z a [r a xs e ][r a x/x Q ] 



-a 



(r a s e )(r a /x a ) - (l/3)r2(s e /0 



2m 2 2m 2 
Among those terms, only V6 and Vs contribute to &f- 



2 (ge/Q 

3 d 



(13a) 
(13b) 



(14a) 
(14b) 



(14c) 
(14d) 



C. Contact type NRQED interactions. 



Here we introduce corrections already mentioned in the introduction, which enter into the NRQED Lagrangian as 
contact type interactions, since they reproduce effects of the relativistic scale. 
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Radiative interactions of order a(Za)Ep: 
8irZ fi p f 13 3 



and of order a(Za) 2 \n 2 (Za)Ep 



• Zemach term ((Za)(m/A)Ep): 

V Z = -2aR z {\ + 55*) |^|| (s e • I a )S(r a ). (17) 

where R z = 1.045(16) @]. 

• Recoil correction of order (Za)(m/M)Ep [(j [Io|: 

V recoil = a [5.48(6) • 1(T 6 ] i^i^- (s e ■ I a ) J(r a ). (18) 

The correction of order (Za) 2 (m/M)Ep from Eq. (J3J) has been omitted since it has to be considered in the 
context of higher order corrections, which are out of the scope of present consideration. 

• proton polarizability (llj : 

V pol = a [1.4(6) • 10- 6 ] JjL|p. (s e • I a )8(r a ). (19) 

D. Perturbation formalism. 

To calculate the bound state problem we use the nonrelativistic Rayleigh-Schrodinger perturbation theory, where 
the starting point, the zero order approximation, is the nonrelativistic Schrodinger equation: 

Hq^o = £ *o, (2°) 
and the perturbation is the effective Hamiltonian H e g derived from the NRQED Lagrangian and 

d 



AE = (H oS ) + (H cS Q(E - Ho^QHcs) + (H cS ) 



dE HcS 



+ ... (21) 



where Q is a projector operator on the subspace orthogonal to the zero-order wave function. H e g has contributions 
of different orders in a: 

H cS = + + + ... 

Then for our case the complete contribution at order of a 6 (m/M) to the hyperfine structure of hydrogen atom and 
ion can be expressed by 

AE Ms = (H^) + (h^Q(E - Hor'QH^) . (22) 

where and are parts of the Breit-Pauli Hamiltonian taken so that the second term in (|22|) contributes to 
that particular order. Since the effective Hamiltonian for the HFS does not depend on E explicitly, the last term 
of Eq. (j2Tj) vanishes. In the following, the second-order contribution and the first-order contribution (H^) will be 
denoted AEa and AEb, respectively. 
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III. HFS IN THE HYDROGEN GROUND STATE 



In the remaining part of this work we will be using the atomic units (e = % = 1 and c = a -1 ). First, we consider 
the case of the HFS of the ground state of a hydrogen atom. Our derivation is somewhat similar to the one done by 
Nio and Kinoshita in [l8l |. The divergent part is, however, treated in a different way by explicitly separating out and 
cancelling the divergences. We start from the nonrelativistic Schrodinger equation: 



where 



V 



Z 
r 



(23) 



(24) 



A. Separating divergences in the ma 6 (m/M) order effective Hamiltonian 

The effective Hamiltonian of order ma 6 (m/M) is obtained from Eqs. (fT2| and (fT4| . interactions V4 and Vs, expressed 
in atomic units. It has the form: 



tfW = a 4 



3m e 



Z I 



Am 2 



2m e r 4 



where \x p = (iip/M p ) l p , then using the relation 



(Z 2 /r A ) 



(PeVY 



pIv 2 + v 2 P i 



+ (peV 2 p e ) - AnZ (VS(r)) 



(25) 



(26) 



obtained by integration by parts, and equation p 2 ^o = 2(Eq — V)^o, one gets: 



AE { V = (H^) = a 4 



Hp 



3Zm e M p 



2ttZ(V6(v)) + . v „ _ £ /y2\ + (PeV 2 Pe) inZEo (8(r)) 



2m P 



(Se-lp). (27) 



In Eq. (|27|) . the divergent contributions are now explicitly collected in the first two terms. 

B. Separating divergences in the second order contribution 

The second order contribution of order ma 6 (m/ M) to the spin-spin interaction can be easily identified from various 
combinations of terms of the Breit-Pauli Hamiltonian and may be written: 



AEf = 2a 4 



Z 



8m5* 8m 2 



47r<5(r) 



Q(E Q - H )- 4 Q 



3m. 



(28) 



This contribution is divergent due to presence of the delta-function operators on both sides of the second order 
iteration. 

Let us consider the two operators: 



H 



(i) 



47r5(r), 



H 



(2) 



8m? 8m?. 



Air6(r), 



and introduce the wavefunction solution of equation 



(E - H )*$ = Q [4rr*(r)] * = QH%>* . 
behaves as 1/r at r — > 0. We introduce a less singular function defined by 



(i) 



2m e * 



C/ 1 * + ^ 1) , 



(29) 



(30) 



where U\ = —^r- — ^r-V . The function ^ B ' behaves as lnr at r — > 0. It satisfies equation 

(e -h )¥^ = (h^-{h^))^, 



where 



H' B {1) =-(E - H )U X - U X {E Q - H ) + H$ . 



(31) 



(32) 



Similar computations can be applied to the scalar part of the Breit-Pauli Hamiltonian H B 

- 4 Ztt 



(2). 



(E - H )^ = Q 



P 



8njp 2ml 



(2) _ gjW , ,T,(2) _ 



4m„r 



S(r) 



(2) 



Uo = 



z 



1 



(33) 



4m e r 



4m„ 



-V. 



and 



H 



'(2) 



-(So - H )U 2 - U 2 (E - H ) + H% 



(2) 



(34) 



Using systematically that 



/ 



H B Q(Eo — Ho) QH B 
one may separate the divergent singularities in the following way 

.4 4 Mp /„ T V / lT ,_ jjm^tT?. TT_\-lrMI^) 



or < * 



,(2) 



QH B 1] 



3m e M p 

2\i v 
3m e M p 

3m e M p 

2^p 
3m e M p 



s e -I p ) (*o H B Z) Q(E Q - Hoy'QHg 1 



(s P -I 



*> (2 (*0 



(2) 



Ui 



(2) 



(2)\ 
5 / 



(se-I P ) (2(* | (H^ - (h { b ] )) U 1 \^o)+2(^o\h ( b ) Q{Eo-Ho)- 1 QH' e 



(i) 



(35) 



(S e -I p ) 2(*o 



rr(2) 

2/* 



S 2) ) 



2 *r 



C/ 2 (^ (1) - <i^ (1) }) 



H' B 2 ^Q{Eo — Ho) 1 QH' B 1 ^ 



The first two terms of the last expression may be rewritten as the average of a new effective Hamiltonian contributing 
to the ma e (m/ M) order: 



ff'(6) 



t 2/ip 
3m e M p 



(i) 



t/ 2 



-Ui(E - H )U 2 - U 2 (E - H )U! \ (s e -I p ) . 



(36) 



Using regularization and integration by parts in a similar way as in Appendix B of [l3| its expectation value may be 
finally written in the form 



3Zni e M p 



-4m P ( H 



r(2) 
l B 



Z{H B 1] ){V) 



(V) 



2m P 



(37) 



(Se-lp) ■ 



All the divergent terms of Eq. (|35|) are collected as the first two terms of Eq. (|37|) . They clearly cancel out those of 

Eq.ua 

The remaining part of the second order iteration contribution (the last term in Eq. (|35p) is finite 



H' B {2) Q{Eo - Ho)- 1 QH' B {1> *„) (s e -I p ) , 



(i) 



(38) 
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where H'J- 1 ^ and H' B 2 ^ are denned above in Eqs. (|32p and respectively. 

Summing up (H^) and (ff'( 6 )) from Eqs. (|27|) and (|3T[) one gets a finite expression as well 



2^ p 



3Zm e M p 



+ 2E Q (V 2 ) 2E 2 (V) - ^ 4. (S(r)) 



m e 



Am e (H {2) )(V) 



Z(hW)(V) 



2m e 



(39) 



(S e -Ip) 



C. Calculation of expectation values and final result 

We now check that Eq. (|3"9")) leads to the usual result in the case of an f s hydrogen atom. Here, = 2Z 3 / 2 e~ Zr is 
the ground state wave function. First we look for a solution of equation: 



(E - H )¥£ = (H' B (lJ -(H' B (± '))*o, H' B ^> = ~d r - 2U 1 (E -H I 



(i) 



and get 

#W = 4Z(\nr - 1 + Zr)^ . 
The next step is to calculate the expectation value of 



H' B m Q 



(i) 



h' {2) -(h' {2 



(i) 

B /> 



and to get the finite part of the second order contribution 



H' b 2 ^Q{Eq — Hq) 1 QH' B 1 ^ 



3(Za) 



-E F (s e -Ip) . 



The expectation values of the operators involved in AE B are (pl^ 2 p) = Z 2 {V 2 ), (V 2 ) = 2Z 2 , and (V) = —Z. 
This contribution is immediately obtained to be 0, indeed 



AE 



,( 6 ) _ 2Zm e a 4 n p 
3 M n 



2Z 5 -2Z 5 + — + 2Z 5 - — -2Z t 



(Se-Ip) =0. 



Thus, the total contribution is 

AE^ = AEf ] + AEf ] = \{Za) 2 E F <s e -I p ) . 
that exactly matches the well-known Breit relativistic correction [l5[ . 



IV. HYDROGEN MOLECULAR ION H+ 

Now we are ready to study the hydrogen molecular ion. As in the previous section we start from the nonrelativistic 
equation with the Hamiltonian: 



Hn = 



P 

2m e 



V, V = 



Z\ Z 2 



(40) 



We will assume here that Z\ = Z2 = Z and fi 1 = (fj, p /M p )li, fx 2 = (Hp/Mp)l2, where Ii and 1 2 are the two proton spin 
operators. 

The second order contribution of the spin-spin interaction of order ma 6 (mj M) is expressed by 



AE A = 2aH- P 



Z 



8m^ 8m 2 



4tt [(S(n) + (J(r 2 )] 



Q{E - H^Q 



Stt 



3m c 



s e [n 1 8(ri)+fj, 2 S(r 2 ) 



(41) 
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The effective Hamiltonian of order ma 6 (m/M) is obtained from Eqs. (|T2|) and (fT4|) of Sec. II. Now we have three 
interactions, V4, Vq, and Vs, because we have as well the seagull interaction Vq with two different nuclei, 



V 4Q = 



1 



4m2 



Pe- 



3 r° 



4 Z f 2(nr 2 )(s e /ij) (rir 2 )(s e /i 7 )-3(ris e )(r 2 /x 2 )-3(r 2 s e )(ri/x 1 ) 



6m; 



r, r: 



V 8a = a 4 



6m; 



2 ( S eMa) , ^(s e /X )-3(r S e )(r Q /i Q ) 



(42a) 
(42b) 
(42c) 



where fij = fj, 1 + /x 2 . It is convenient to separate the effective Hamiltonian into two terms: scalar and tensor, 



Hp 



3m e m p 



4 ^{rf,47r[5(r 1 )+5(r 1 )]} + 



Z / 1 1 2rir 2 



2m e ^r* r 2 r^r? 



(Se-I) 



6m« 



r i( s eM 2 )-3(riS e )(ri/x 1 ) ^(s e /x 2 )-3(r 2 s e )(r 2 /x 2 ) 



(rir 2 )(s e /x j )-3(ris e )(r 2 /x 2 )-3(r 2 s e )(ri/x 2 ) 



(43) 



h[ has a finite expectation value, and since it does not contribute to bp , its consideration will be omitted in what 
follows. The divergent terms are encountered only in the scalar Hamiltonian. 



A. Separating divergences in the second order contribution 



The operators which appear on the left and on the right of the second order iteration are: 

^ 4 Z 



H 



(1) 



Now we have to separate the singular part using the method outlined in the hydrogen case. We set : 



H 



(2) 



(E -H )^ B ' = Q 47r(<S(n) + a(r 2 ) 







¥ B } = 2m, 



H'^ = -(E - H )U! - U^Eo - H ) + H B A 



1 1 

n r 2 



(i) 



Similarly, one gets for H B 



(2). 



(E Q - H )^ = Q 



^ = 



Z 



1 1 

ri r 2 



Ui = 



1 



4m e 



-V. 



Hf ] = -(E - H )U 2 - U 2 {E Q - H ) + Hf 
Applying these transformations to the second order iteration term we arrive at 



AE A = a 4 — ^ 

3m e M p 



^ 2) -<4 2) >)^ 



U 2 {H' B {1) - (H' B {1) )) 



H'g 2 Q(Eq — H ) 1 QH' B 



Se-I) 



Again we pick out the first two terms which can be recast in the form of an effective Hamiltonian: 



3Zm e M p 



p 4 V+Vp 4 (Vp 2 V) , , / (2)\ Z(H ( n ] )V 

1 y y * -V+E V— 4m e (H l B ' )V J X B 1 



477j2 



2m P 



2m„ 



(s P -I) 



(44) 



(45) 



(46) 



(47) 
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R (in a.u.) 



FIG. 2: Adiabatic effective potentiai for the reiativistic correction of order ma 6 (m/M) to the spin-spin interaction coefficient 
6f- Energy scaie for A6f is in (a 4 n P /M p ) x (1 a.u.). 









L = 1 


L 


= 3 


w 


new 


m 


new 


V 


= 


922.9918 


922.9168 


917.5911 


917.5167 


V 


= 1 


898.8091 


898.7371 


893.7545 


893.6831 


V 


= 2 


876.4542 


876.3851 


871.7277 


871.6592 


V 


= 3 


855.8124 


855.7460 


851.3984 


851.3325 


V 


= 4 


836.7835 


836.7197 


832.6682 


832.6049 



TABLE II: Results of numerical calculations for low ro- vibrational states of the spin-spin interaction coefficient &f (in MHz). 



Its expectation value can be rewritten as follows: 



3Zm e M p 



_ {v * } + <p|£p> + 3Eo (v * } _ 2E , (v) 



-4m P ( H 



r(2) 
l B 



Z(H™){V) 



(V) 



2m P 



(s R -I) 



and the divergent terms, the first two terms, are now written explicitly. 



(48) 



B. Removing divergences and final expressions 

The remaining part of the second order iteration contribution (last term in Eq. (|46p is now finite 



3m e M p 



H' B (2) Q(E - H y l QH' B {1) ¥„) (Be- 1) , 



.« and H' B {2) 

Summing up \H^\ and (iJ'( 6 )) we get a finite expression as well 



where H B ' and H' B ' are defined above 



B 3Zm e Mp 



+ 2E Q (V 2 ) 2E 2 (V) ^ 4nZ (£( ri )+£(r 2 )) 



-4m e (h% ) )(V) + 



Z{H B 1] ){V) 



2m e 



(49) 



(50) 
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V 


n l 


this work 


experiment 


1 

4 


ooO. f o4 


ood. / z\y 


ooD. / zy 


5 


819.280 


819.219 


819.227 


6 


803.227 


803.167 


803.175 


7 


788.558 


788.501 


788.508 


8 


775.221 


775.166 


775.172 



TABLE III: Comparison of the spin-spin interaction coefficient bp (in MHz) with experiment. Ref. pi is the Breit-Pauli 
approximation with account of the electron anomalous magnetic moment. L = 1. 



contribution v = 4 « = 5 





836.7835 


819.2801 


(Z«) 2 


0.0510 


0.0511 


a(Za) 
a(Za) 2 ln 2 (^a) 


-0.0804 


-0.0787 


-0.0067 


-0.0065 


ASz 


-0.0335(5) 


-0.0328(5) 




0.0049(1) 


0.0045(1) 


A_Bp i 


0.0012(5) 


0.0011(5) 



b F (new) 836.7197(10) 819.2187(10) 



TABLE IV: A summary of contributions to the spin-spin interaction coefficient 6f (in MHz). 

V. RESULTS AND CONCLUSION 

Results of numerical calculation as a function of a bond length for the relativistic correction to the HFS within the 
framework of the two-center problem are shown on Fig. [2] In our study we use the variational exponential expansion 
introduced in [lj|. In fact, the adiabatic effective potentials from Eq. ([50)) have already been obtained in the previous 
work [l3| and only the second order perturbation term (Eq. ([49)) ) with modified operators H'^p and H'^ require 

some additional numerical efforts. The potential of the total effective Hamiltonian AE'g^ tends to zero when R — » 0, 
or R —* oo, as it may be expected from the analysis of the hydrogen atom ground state HFS. 

The relative numerical accuracy of the potential curve plotted on Fig. [5] is estimated to be ~ 10~ 5 , however the 
adiabatic approximation itself limits the final uncertainty of the relativistic contribution of the ma e (m/M) order to 
the spin-spin interaction coefficient bp to be about 0.1 kHz (3-4 significant digits in Abp). The other contributions 
which are described by Eqs. (JTSJ)— (jTSJ) may be obtained using the previously calculated mean values of the delta 
function operators [l^]. The final results for the new theoretical value of the coefficient bp for the low ro- vibrational 
states are presented in Table [TTJ 

An experimental value for bp can be uniquely calculated by using the mixing parameters (20| of the states (F, J) : 
(1/2, 1/2) «-> (3/2, 1/2) and (1/2,3/2) <-> (3/2,3/2), to restore the structure of pure F = 1/2 and F = 3/2 multiplets 
and then take a difference between statistically averaged splittings of these multiplets. In Table Hill a comparison with 
experiment is given. As it may be seen the newly obtained results improve the agreement with the experiment by 
about a factor of 6. The error bars for bp from the experimental data are to be about 1 kHz as it follows from the 
claimed accuracy of Ref. Q. On the other hand from the comparison with the hydrogen atom case, the theoretical 
uncertainty should be no more than 2-3 kHz. That indicates substantial discrepancy between theory and experiment 
of about 6-9 kHz. 

In order to try to explain this discrepancy we have checked several effects which may have impact on the spin-spin 
interaction. The leading order retardation effects in the nonrelativistic interaction region (2lj as well as the cross 
terms of the second order perturbation (when electron interacts with both protons in Hj) for the proton structure 
dependent contributions are estimated cither equal to zero or negligibly small. We have also analyzed the effect of 
the g/u symmetry breaking which is essential for high v states, say, for v = 19, it leads to a few MHz shift in energy 
[22| . However for the states below v — 10 this effect is smaller than 1 kHz and, thus, the gap between theory and 
experiment can not be accounted for by the g/u mixing. A possible explanation is that the higher order corrections 
(a (m/M)Ep and a 3 Ep) may give significant contribution. 

In conclusion, the consideration of ma 6 and (partially) ma 7 order corrections, as well as proton finite size effects, 
has allowed to improve significantly the agreement with experiment, to about 10 ppm. The remaining discrepancy is 
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somewhat larger than expected from comparison with the hydrogen atom, and further theoretical work to improve 
the HFS intervals is needed. In any case, a new independent experiment is highly desirable. 
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